ON TRAJECTORIES OF VORTICES IN THE COMPRESSIBLE 
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Abstract. For the model of a compressible barotropic fluid on a two dimensional rotating 
Riemmanian manifold we discuss a special class of smooth solutions having a form of a steady 
non-singular vortex moving with a bearing field. The model can be obtained from the system of 
primitive equations governing the motion of air over the Earth surface after averaging over the 
height and therefore the solution obtained can be interpreted as a tropical cyclone which is known 
as a long time existing stable vortex. We consider approximations of I- plane and /3 - plane used in 
geophysics for modeling of middle scale processes and equations on the whole sphere as well. We 
show that the solutions of the mentioned form satisfy the equations of the model either exactly or 
with a discrepancy which is small in a neighborhood of the trajectory of the center of vortex. We 
perform a numeric study of the change of the shape of the vortex affected by the neglecting the 
discrepancy term. 

Key words, mathematical model of atmosphere, exact solution, vortex trajectory, tropical 
cyclone 
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1. Introduction. A lot of physical models can be reduced to the equations of 
motion of compressible medium on a two-dimensional manifold. One of the most 
important is the model of dynamics of the atmosphere. The vertical scale of the 
atmosphere is much more larger than the horizontal one, therefore one can reduce 
the initial primitive three-dimensional system of equations to a two-dimensional one, 
convenient for describing many kinds of motions of middle and global scale, using a 
special procedure of averaging over the height. 

There are a lot of papers where the equations governing the atmospheric motion 
are reduced to the two-dimensional system of incompressible Navier-Stokes equations 
with viscosity (see [31], [12], and references therein), on this way a lot of results on 
existence and uniqueness of solution to the Cauchy problem were obtained. The effect 
of compressibility makes the situation more difficult and the results on the global in 
time correctness for the viscous case are very scared (a recent state of art can be 
found in [14). If we neglect the viscosity as they usually do in the meteorology (thus 
we get the quasilinear hyperbolic system for compressible fluid), we find ourself in a 
more complicated situation since the solution to the Cauchy problem basically loses 
its initial smoothness with time. Recent developments in the study of the Cauchy 
problem for the equations for compressible fluids are reviewed in [S]. 

Nevertheless, it is well known that there exist very stable vortices in the at- 
mosphere that move without changing their shape for days and even weeks. Most 
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notorious example is the tropical typhoon, one of the most devastating weather phe- 
nomena in the world. In practice, the most important problem is to describe (or 
better predict) the path of typhoon, whereas the processes inside the vortex are of 
theoretical interest. 

In our earlier papers [29], [30] we developed a theory on possible trajectories of 
a stable vortex with a linear profile of velocity governed by the system of equations 
of compressible fluid on a plane (see Sec l4.1l for details). These vortices can be inter- 
preted as tropical cyclones in the conservative phase of its dynamics, thought they 
have a number of nonrealistic features. First, the components of the velocity and 
pressure rise unboundedly as the distance from the center increases (we call these 
vortices non-localized). Second, the curvature of the Earth surface was not taken into 
account even in a simplest form. Nevertheless, as it was shown in |30j by comparing 
with observational data, the theoretical trajectories of vortices can imitate the real 
tropical cyclone paths. In fact, they are a superposition of two circular motions and 
correspond to the natural circular (or rather parabolic, taking into account the change 
of the Coriolis parameter with the latitude) trajectories, loops, reversal points, etc. 
Moreover, the explicit expression for the trajectory was obtained. 

In the present paper we show that the trajectories of localized vortices that model 
the physical process more adequately are very similar to the trajectories of the non- 
localized vortex with linear profile of velocity considered in [30] under some realistic 
assumptions on the background field of pressure. 

It is worth mentioning that the exact form of equations describing the atmosphere 
is not known. More exactly, in the primitive system of equations we can take into 
account some addition forces or sources of energy, e.g. due to phase transitions, and 
the exact form of these terms can hardly be found, particulary after the averaging 
procedure. 

We will base on the assumption on existence of special form of solution for the 
system of equations of the atmosphere dynamics (steady vortex moving in a bearing 
field) together with the assumption on a smallness of the source term that guaran- 
tees the existence of this special solution, at least in a domain of space that we are 
interesting to control. 

We can also look at the problem from another side and call a vector-function 
the "<5 - approximate" solution to a certain system in a domain D G M 2 x K+ if the 
discrepancy term Q arising after substitution of this function to the system is less 
than S in the uniform norm. Thus, we can consider the term Q either as some source 
term that guarantees the existence of the exact vortex solution or as a discrepancy, in 
the latter case the vortex solution can be considered as an approximative one. 

The paper is organized as follows. In Section [5] we recall the derivation of a 
bidimensional model of the atmosphere dynamics based on the primitive equations 
for compressible viscous heat-conductive gas in the physical 3D space. In Section |3] 
we develop a general approach to the solution in the form of " frozen" vortex moving 
in the exterior field. In Section|4]we demonstrate this method for the vortex solutions 
in the I - plane and the j3 - plane models and consider different shape of possible 
vortices (both localized and non-localized). Moreover, we compare the trajectories 
obtained here with the the trajectory of vortex with linear velocity profile where the 
exact analytical result is available. Further, we perform direct numerical computation 
of the moving vortex in the 2D models of the Z-plane and /3-plane and compare the 
position of the center of vortex with the result obtained by our method. In Section [3] 
we consider the case of a sphere. We construct global exact solutions for non-rotating 
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sphere and study the trajectories of approximative vortices for the rotating case as 
well. Finally, a discussion about our topics is provided in Section [5J 

2. Bidimensional models of the atmosphere dynamics. We do not dwell 
at the procedure of averaging and refer to our paper [30] , where we use the approach 
by [51] and [2] ■ The initial (primitive) three-dimensional system relates to the motion 
of compressible rotating, viscous, heat-conductive, Newtonian polytropic gas [T7] , 
the resulting two-dimensional system consists of three equations for density g(t,x), 
velocity U(t, x) and pressure P(t,x) (the equation for the velocity is vectorial): 

g(d t U + (U • V)U + ILTJ) + VP = Pi, (2.1) 
d t Q + V-( g U)=0, (2.2) 
d t P+(U,VP)+ 7 PdivU = P 2 , (2.3) 

where Pi and F 2 are some source terms (see [30] ) for details, L = y ^ J J , i if 

the Coriolis parameter, a smooth function of x, its concrete form we discuss below. 
In fact, system (|2.1[) - (|2.3p consists of equations of balance of mass, momentum and 
total energy, t £ R + , x € A4, where M. is a 2D Riemmanian manifold. 

We can consider the system at any 2D Riemannian manifold A4 (with or without 
boundary). Thus, we have to take into account the curvilinear metrics of the space 
and write the system in the respective curvilinear coordinates. In this case we mean 
by partial derivatives the covariant derivatives |37j . In particular, to model the motion 
of planetary scale we should use the spherical coordinates ([25 ). Nevertheless, for the 
phenomena of middle scale it is convenient to use the approximation of plane, without 
taking into account the effects of curvature. 

Now we recall shortly the procedure of averaging. Let p, u — (1*1,112,1*3), p, 

denote in the three-dimensional density, velocity and pressure. Namely, all these 

functions depend on (t, x, z), z 6 M + . Let us introduce <f> and / to represent for 

taking the average of <f> and / over the height, respectively. The averaged values are 

/•oo y r°° 

introduced as follows: (f> :— / (f>dz, f := — / pf dz, where <fi and / are arbitrary 

Jo P Jo 

functions, and denote g(t, x) = p, P{t,x) = p, XJ(t,x) = (1*1,1*2)- Moreover, the usual 

adiabatic exponent, 7, is related to the "two-dimensional" adiabatic exponent 7 as 

27 - 1 
follows: 7 = < 7. 

1 

The impenetrability conditions are included in the model. These conditions ensure 
that the derivatives of the velocity equal to zero on the Earth surface and a sufficiently 
rapid decay for all thermodynamic quantities as the vertical coordinate z approaches 
to infinity. In other words, the impenetrability conditions make sure the boundedness 
of the mass, energy, and momentum in the air column. They also provide the necessary 
conditions for the convergence of integrals. 

The source terms Pi and P2 include basically the (turbulent) viscosity, heat con- 
ductivity, and may be something else. In fact, the form of Pi and P 2 for real geophys- 
ical model depends on the model considered, in any case these terms are believed to 
be small and to have a regularizing effect. Indeed, the solution to system (|2.1[) - (|2.3p 
(i.e. the compressible Euler system) generally loses smoothness within a finite time 
with formation of shocks, and it needs to introduce some smoothing terms to avoid 
this phenomenon not inherent to the atmospheric motion. 
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For three-dimensional equations it is the practice to introduce the entropy func- 
tion s, which is connected with the pressure and density through the state equation 

p = a e s p 1 , 

with a constant a > 0. 

We can also introduce the 2D entropy S(t, x) connecting with g and P similarly 
to ([2]), where we use 7 instead of 7. Then instead of equation (|2.3[) we get the equation 
for the entropy 

dtS+{U.yS) = F 3 , (2.4) 

with a source F3, see [30] for details. 

For the sake of simplicity and to make our idea clear in this paper we will consider 
only the barotropic case of constant entropy So, in other words, the state equation is 

P = Cpi, 7 > 1, C = const > 0, (2.5) 

and (|2.4p is assumed to hold identically with F3 = 0. 

Thus, the system under consideration is reduced to two equations (|2.1[) . (|2 . 2[) . 

1-1 

It will be convenient to introduce a new variable tt = P -< . For the new unknown 
variables 7r(t, x), U(i, x) we obtain the system 

d t U + (U • V)U + ILV + c Vtt = F, 
d t -K + (Vtt • U) + (7 - 1)tt divU = 0, 

with a new source term F(t,x) and Co = ^j-C~. 

3. "Frozen" vortex: a general approach. Let us change the coordinate sys- 
tem in such a way that the origin of the new system x = (27, X2) is located at a point 
X(f) = (Xi(t),X2(t)) (here and below we use the bold font for x to denote the local 
coordinate system). Now U = u + V, where V(t) = (Vi(t),V 2 (t)) = (X^t), X 2 (t)). 
Thus, we obtain a new system 

d t u + (u • V)u + V + |(X(t) + x) L (u + V) + c = F, (3.1) 



d t n + Vtt • u+ (7 - l)7rdivu = 0. (3.2) 

Given a vector V, the trajectory can be found by integrating the system 

X 1 (t) = V 1 (t), X 2 (t) = V 2 (t). (3.3) 

Below we look for a solution with special properties, namely, a steady divergency 
free vortex " frozen" into a certain exterior pressure field. Together with this we study 
a class of sources F allowing the existence of such solutions. Let us stress that we do 
not prescribe specific boundary conditions and only assume that 

u does not depend of t, u(0) = and div u = 0. 

This means that there exists a potential $(iri, x 2 ) such that 



(3.4) 
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Thus, the velocity held satisfies the Cauchy-Riemann conditions 

d xi ui = d X2 u 2 , d X2 ui = -d Xl u 2 . 

In [TT] it is proved that for the case of Z-plane this condition is necessary for the 
existence of a stable vortex. 
Further, we assume that 



where 



7T = T7 (X 1 ,X 2 ) + 7T 1 (t,X 1 ,X 2 ), 



(V^ ,Vi$) =0. 



(3.5) 



(3.6) 



We call the time-dependent component m a bearing field of pressure. It will be 
important for us that the gradient of the bearing field is rather small. The time- 
independent part of pressure ttq relates to the vortex itself. 

Let A be an appropriate (2 x 2) matrix which depends on the situation to be 
discussed below. We consider the following equation for the potential of velocity $ 
and the steady part of pressure ttq: 



(V_l$, V)V_L$ + AV_L$ + CqVtTo = 0. 



(3.7) 



Let us suppose that we succeed to solve system (13.71) . (|3.6[) . Then from (I3.2p we 
get a linear equation for -k\. 



(3.8) 



V(t) + l{X(t) + x)LV(t) + (l(X(t) + x)L - A)u(x) + c V^i(t, x) = F, (3.9) 



d t TTl + V7Tl • U = 0, 

which can be solved for any initial condition 7Ti(0, x) 
Further, from (13.11) we obtain 



or 



X(t) + l{X{t) + x)LX(t) + (l(X(t) + x)£ - A)u(x) + coVvr^i, x) = F. 
Let us denote 

Q = -(l(X(t)) - l(X(t) + x))iX(t)+ 

V7Ti(t,x) . 



(l(X(t) + x) L - A)u(x) - c V7ri(i,x) 
From f3~TUl) and (|3TTT|) we get 

X(t) + l(X(t))LX(t) + coVwx(t, x 



x=0 



(3.10) 



(3-11) 



= F-Q. 



(3.12) 



If F equals Q, then the position of the center of vortex under consideration can be 
found from the following equation: 



X(*) + ;(X(t))XX(t) + coV7ri(t,x) 



0. 



x=0 



(3.13) 
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Thus, (|3.1|) can be rewritten as 

(u- V)u + J 4u + c V7r + 
X(t) + l(X(t))LX(t) + cV^x) = (3.14) 

F-Q, 

the terms in the first row depend only on the space variable x, the terms in the 
second row depend only on t, therefore for the case F — Q = we obtain a complete 
separation of variables for solution satisfying (|2.6p . (13.41) . (|3.5[) . (|3.6|) . (|3 . 8[) . (|3.13|) . 

If F — Q is not zero, however it is small (for example, the uniform norm \\F — Q\\ < 
5), we can talk about a "5- approximate" separation of variables, F — Q plays a role 
of the discrepancy. 

Definition 3.1. We call a couple (n(t, x), U(£, x)) the 5 - approximate solution to 
system (|2.6p in a domain D{8) £ M+ xAi, if(w,\J) satisfies (|2.6p with a discrepancy, 
whose uniform norm is smaller than 5 for any (t, x) £ F>(S). 

In fact we have to choose the matrix A in such a way as to ensure the smallness 
of the discrepancy term F — Q. Since we are interested in studying the behavior of 
solution near the origin of the moving coordinates system, we expand in (|3.1ip the 
known functions Z(x), u(x) and itiit, x) in the origin and get Q = 0(|x|) provided 
X(i) keeps boundedness. Basically a solution to the nonlinear equation (|3. 12|) blows 
up in a finite time i*, this entails an unboundness of X(i). Thus, our considerations 
are valid for t < , until the formation of the shock wave in a neighborhood of the 
trajectory X(i) of the moving coordinate system. 

Let us list once more the steps of our method. First we should solve (I3.6p . (|3.7p , 
then set an initial distribution for the bearing pressure field 7Ti(0, x\, X2), then find 
the respective time-dependent pressure component 7Ti (£, x\ , X2 ) and at last integrate 
(|3.10p to find exact (for F = Q) or approximate (for small F — Q) trajectory of the 
exact (or S- approximate) vortex solution. 

Thus, the main problem is to solve the system (I3.6p . (I3.7p . with respect to un- 
known scalar functions $ and ttq (may be under suitable boundary conditions). 

Let us find a necessary condition for the function $ that allows it to be a part 
of solution of (|3. 61) . p. 71) . We take the inner product of (|3.7p and and get the 

master equation 

V Ki:C2 $((V ;ci $) 2 -(V X2 $) 2 ) + (AV ± $,Vi$) = V^SV^SCVx^S-Vxx,!*). 

(3.15) 

If the solution of (|3.15p satisfies the identity 

Vx((Vi$,V)Vi$ + AVi$)=0, (3.16) 

then one can find ttq such that the couple ($, tto) solves Q3.6p . (|3.7p . 
Let us summarize the results of this section. 

Theorem 3.2. Let $(x) be a solution to equation fl3.15p . satisfying condition 
(|3.6p with a matrix A with smooth coefficients dependent only on space variables. 
Further, let 7To(x) be a solution to (|3.7p . 7Ti(t,x) be a solution to (|3.8p with u = Vj_$, 
u(0) = 0, and X(t) be a solution to (|3.13l) . Then the couple (7T,U), where 

n(t,x) =7r (x-X(i))+7r 1 (i ) o;-X(i)) ) U(f, x) = X(t) + u(x - X(t)) (3.17) 

sofoe system ()2.6[) wii/i F = Q, Q is given by (|3.1ip . 

Theorem 3.3. We assume F = 0. Lei $, 7r , 7ri, X(t) 6e as m Theorem WB. 
moreover, they are classical solutions to the respective equations. Assume |X(f)| io 6e 
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bounded for t < T . Then there exists a neighborhood D(S) of the trajectory x = X(t) 
such that (|3.17p is the S- approximate solution to system (12.61) for t < T . The dimen- 
sions of D(S) depend on derivatives of known functions, namely, |VZ(x)|, |V7r(t, x)\, 

2 

|h|, where h = (hi,h 2 ), hj — J2 (IL - A)ijUi, j = 1,2, u = (ui, u 2 ). 

i=l 

Proof. The proof of Theorem 13.31 follows from estimating the term Q given by 

Below we dwell on important for geophysics cases of application of the method. 

4. Steady vortices on a plane. The "plane" models are used in meteorology 
to describe the processes of small and medium scale, where the curvature of the Earth 
surface does not play a crucial role and the change of the Coriolis parameter can be 
modeled in a simplest way. 

Let Xo = (£01,2:02) be a point on the Earth surface, ipo be the latitude of some 
fixed point Xo. The simplest possible "plane" model is the so called Z-plane model, 
where the Coriolis parameter / is treated as a constant: I = Iq = 2Qsin<poj f2 is the 
vertical component of the angular velocity of the Earth rotation. The model, that is 
believed to be more adequate to describe the weather processes, is the /3-plane model, 
where the Coriolis parameter I is approximated by a linear function: I = Iq + 13x2, the 
constant j3 = ^-cosfo, where R is the radius of the Earth. In fact, in the Z-plane 
model we neglect the parameter /3, which is much more smaller than Iq. For example, 
for ip a = 30°, Z « 7.3 x 10" 5 s _1 , /3 « 2 x 10 -n s -1 , where we used the values 
R = 6.4 x 10 6 m and uj 3 = 7.3 x lCT 5 rad/s. 

It can be readily checked that for I = Iq = const every potential 

$ = $( r ), r = ^Jx\+x 2 2 , 

solves the master equation (|3.15[) with A = IqL. We can get different shapes of vortices 
choosing different solutions. 

For the /3-plane model we will use the master equation with the same matrix 
A = IqL. Recall that the choice of this matrix can be different, nevertheless it is 
caused by a desire to make the discrepancy term F — Q as smaller as possible. 

4.1. Example 1: steady vortex on the I— plane (linear profile of veloc- 
ity). We begin from the simplest case where for F = Q = we have the complete 
separation of variables in equation (|3.14[) . Here we can obtain an exact solution. In- 
deed, if we choose $ = — {x\ + x\), bo = const, we get the velocity field with a linear 
profile 

m = box 2 , u 2 = -b xi. (4.1) 

Further, if we choose the initial data for 7Ti as a linear function of the space variables, 
7ri(0,x) = Mio^i + M 2 qx 2 + K Q , A/10, M 2 q, K = const, we get 

vi(t,x) = Mi(t)sx + M 2 (t)x 2 + K(t), (4.2) 

where 

Mx{t) = M10 cos&o* + M 20 sin6 t, M 2 (t) = M 20 cosb Q t - M w sinb t. (4.3) 
The trajectory obeys the equation 



X(t) + i iX(t) + coM(i) = 0, 



(4.4) 
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where M(t) = (M^t), M 2 (t)). 

The solution can be found explicitly: 

X l{ t)= X l( 0) + « + ^ 

L UQL 

fVi(0) c Q M 20 \ . /V a (0) , c M 10 , 

c M 20 . , . , c Mio 

sin oqi + 7—77 jT cos oqi, 



X a (t) = X 2 (0) 



b (bo-l) 6o(6o-0 
Vi(0) , c Af 20 



,V 2 (0) c M 10 \ , /V1(0) c M 20 , 
co^io . , . , CqM 20 

■ Sin Ogl + - — — COS OqI, 



6o(6o-0 u 6o(6o-0 

(for I 7^ &q). This is a superposition of two circular motions: the respective frequencies 
are 7^- and i~. 

On this way we obtain the same solution as in [29 , 30 . The above formulae 
were used in [3D] to imitate and forecast the trajectories of tropical cyclones based on 
real observational data. In fact, we are going to show that a further sophistication of 
model do not implies significant difference in the behavior of the vortex trajectory, at 
least for the "realistic" (relative to meteorological data) values of parameters. 

4.2. Example 2: approximately steady vortex on the /3— plane (linear 

profile of velocity). Let us choose again A = loL, $ = (x\ + x\) , tti(0,x) = 

M w x\ + M 2 qX 2 + Kq, M w , M 2 q, Kq — const. Thus, (14.11) gives the velocity field 
and (|4.2[) gives the background pressure field %% as in Ex.1. 

To obtain the solution of the form of the steady vortex we need to set F = Q. 
where Q (see p.lip ) has the form 

Q = /?x 2 L(X(t) + u(x)) +/3I 2 (t)Iu(x). (4.5) 

To find the trajectory of the steady vortex we have to integrate the equation 

X(t) + (Z + PX 2 (t)) L X(i) + c M(t) = 0. (4.6) 

If F — 0, we can consider Q as a vector of discrepancy, and due to the smallness 
of (3 for real meteorological models the value of Q is essentially small until the blow 
up time of the respective solution to the ODE 



4.3. Example 3: steady localized vortex on the I— plane. As was men- 
tioned before the velocity field (14.11) rises at infinity and is not realistic. Thus it is 
natural to consider the vortex in some sense localized in the space. We will call the 
vortex "localized", if the velocity profile looks like the observational data (e.g. [33], 
see Fig |4.ip . This kind of initial data is used for the majority of numerical compu- 
tations simulating tropical cyclones. It is catching to select the form of the vortex 
basing on the experimental data, for example, in the class of potential functions 

$ W = n ° 2U ' a > 0, a > 0, fc > 0. (4.7) 
(1 + ar A ) K 
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Figure 4.1. The tangential component of (relative) velocity in m ■ s _1 on different stage of 
the typhoon formation (August 30, August 31, September 2). The hurricane Anita, North Atlantic, 
1977 \34^ - The horizontal scale is in kilometers. 

Indeed, the dashed line in Fig |4.2l presents the graph of the modulus of the radial 
component of the corresponding velocity field for k = |, a — 7.8 x 10 6 m 2 /s, a = 
5 x 10~ 10 m~ 2 . The parameters give a correct maximum of velocity and a realistic 
decay of the velocity at infinity. Nevertheless, one can remember that in our model 
we use the averaged over the height velocity and to construct the graph analogous to 
Fig l4.ll for averaged velocity we need to analyze the experimental data for different 
height levels (that we do not have at our disposal). Thus, the profile of velocity for 
averaged data may be different from those that presented in Fig |4.1l At least, the 
value of amplitude a should be much more smaller due to the different direction of 
rotation of the air for lower and upper parts of cyclone [M]. Moreover, if we use the 
pattern (|4.7p . we cannot solve explicitly equation (I3.8[) for the background pressure 
field. Therefore we use a more localized (exponentially decaying) solution to the 
master equation (|3. 15|) with A = IqL: 

$ = -BoeS (*i 2 +^ 2 ), (4.8) 

here Bo is a constant, a is a positive constant. The respective velocity field 

Ul = B ax 2 e-%( Xl2+X22 \ u 2 = -B a xi e~* ^ +x ^) (4.9) 

vanishes at infinity and has the same asymptotics as (|4.ip at the origin, with 6 = B a. 
The solid line in Fig |4.2l shows the modulus of the radial component of velocity for 
B = 3.5 x l(T 6 m 2 /s, a = l(T 9 m- 2 . 

We are going to show that for the values of parameters relevant in the meteorology 
such strict localization does nor lead to a big difference in trajectories of the localized 
and non-localized vortices. 

It can be checked that condition (13. 16)) is satisfied and therefore we can find the 
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Figure 4.2. The tangential component of (relative) velocity in m ■ s 1 for a localized vortex: 
exponential decay (solid line) and power decay (dashed line). The horizontal scale is in kilometers. 



scalar function ttq- Namely, 

7r (xi,»a) = ^ (±Bl<Te-°&+'%> - l B e~i^ 2 A . (4.10) 
Further, equation (|3.8[) can be explicitly solved, its general solution is 

' - 1 < ' ''arcanf- 

X 2 



7Ti = G ( x\ + x 2 , t - — — e 

-DO CT 



with any differcntiable function G. For example, for initial background field 

7ri(0,xi,x 2 ) = Ro + 4>(xi,x 2 ) (MiqXi + M 20 x 2 + K ), (4.11) 



R , Mio, Mio, K = const, fo(xi,x 2 ) = e - a °^ +x ^\ 

we get 

7Tl(i,Xl,X 2 ) = 

-Ro + 4>q{xi,x 2 ) {Kq + M w (xi cos((t)(xi, x 2 ) B at)- x 2 sm(qb(xi,x 2 ) B crt))+ 
M 20 (x 2 cos(4>(xi,x 2 )B at) + xi sin(</>(xi , x 2 ) B at))), 

with <j>{x 1 ,x 2 ) = e~%^+ x ^. 

The trajectory of the vortex can be found by integration from equation (|3.12p : 



X(t) +l L±(t) + C V7Tl(t,x) 



0. 



x=0 



(4.12) 



It can be checked that 



V7Tl(t,x) =(Ml(t),M 2 (t)), 
x=0 



with the same Mi(i) and M 2 (t) as in (|4.3|) . 

Thus, equation (|4.4p that governs the trajectory of the non-localized vortex for 
this specific background pressure field coincides with (I4.12j) . 
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The discrepancy term (for F = 0) here is 



Q = c 



Vvri(t,x) - V7Ti(i,x) 



(4.13) 



it depends only on the properties of the background pressure field 7Ti. From the 
explicit expression for 7Ti , presented above, one can see that Q is bounded and its norm 
vanishes at every point of plane as the parameters of initial slope of the background 
pressure field M\o and M 2 g tend to zero. 

4.4. Example 4: steady localized vortex on the /3— plane. We use in this 
case the same master equation as in Ex.3 and its solution (|4.8p . Moreover, we use the 
same initial background pressure field tti(0, x) as in Ex.3. The equation (|3~12]) for the 
trajectory of the steady vortex is 

X(i) + (I + /3X 2 (t))LX(t) + cqVttx^x) =0, (4.14) 

X— 

this is the same equation as (14.61) . the discrepancy Q (for F = 0) is 
Q = (ix 2 L (X(t) + u(x)) + [iX 2 {t)L\i{*)+ 

C V7Tl(£,x) - V7Tl(i,x) 

X- 

It can be considered as a small one due to the reasons described in Examples 2 and 3. 

4.5. Comparison of trajectories for the I plane and the /3 — plane 
models. We compared the trajectories of vortices corresponding to approximate 
solutions under values of parameters close to the realistic ones. Namely, we set 
O = 7.3 x 10~ 5 s~ 1 , the geographical latitude </>o = 30°, the Earth radius R = 
6.39 x 10 6 m, c = 0.1 (appropriate dimension), B — 3.5 x 10 3 m 2 /s, a = lCT 9 m~ 2 , 
cr = 8 x l(T 13 m- 2 , M = 2x Kr 3 s"\ N a = lfr 3 s-\ X x (0) = X 2 (0) = 0, 
Vi(0) = -lm/s, V 2 (0) = lm/s. 

The result is the following: the difference between the I - plane and the (3 - plane 
models is very small for these parameters up to 2 days period of computation, further 
the difference increases, however the behavior of the both trajectories remains similar 
(see FiggH). 

4.6. Computational experiments with localized vortices on the plane: 
influence of the discrepancy term. As it was shown before, if the discrepancy 
term is zero, the vortex is steady and it moves according to (|3.12p . Nevertheless, let us 
study the influence of the discrepancy on the shape of vortex. It is natural to expect 
that a big discrepancy destroys the vortex, whereas a negligible one only modifies it. 

Figs l4.4H4.6l present the result of computations for the / - plane and the f3 - plane 
models. Here ip = 30°, l « 7.3xKr 5 s-\ /3 « 2xlCT 11 s- 1 , u 10 = Vx(0) = u 20 = 
V 2 (0) = lOm/s. The initial position of vortex is the origin. The parameter Bq, <t, Co 
are as in the Sec l4.5l The localized vortex is given initially as (|4.10l) . 

ui = u w + B ax 2 e~%( Xl2+X22 \ u 2 = u 20 - B a a x x O^^ 2 ) . (4.15) 



Fig l4.4l corresponds to the pressure for 7Ti(0,xi,X2) = (zero bearing field and dis- 
crepancy), the first picture presents the level lines for the initial data (common for 
all figures), the second one corresponds to the vortex computed for approximately 4 
days. Fig l4.5l corresponds to the vortex on the I - plane in the bearing field. The 
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Figure 4.3. Trajectories of vortices on l-plane (upper curve) and - plane for 7 days period. 
The space scale is in kilometers. 
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Figure 4.4. The l-plane case without bearing field, initial vortex, T = and T 96h. 



parameters of "inclination" of the bearing field are Mio = — 10~ 5 , A/20 = 10~ 5 (weak 
bearing field) for the first picture and Mio = — 10~ 5 , M20 = 10~ 5 (stronger bearing 
field) for the second picture. Computations are made for one day, the space scale is 
the same for both pictures in Fig l4.5l and for the second picture in Fig. 14.41 We can see 
that the vortex preserves its initial shape for the null bearing field, it modifies its form 
and enlarges in a weak bearing field and tends to disappear in a strong bearing field. 
Fig. 14.61 corresponds to the initial bearing field of form (|4. 1 1[) with <tq = 10~ 13 (here 
and below the dimension is appropriate) for the I - plane and the /3- plane models. 
The parameters of "inclination" of the bearing field are Mio = — 10 -5 , M20 = 10 -5 
(weak bearing field, approximately 2 days period). The blue asterisk corresponds to 
the theoretical position of vortex computed according to (I4.12j) and (|4.14[) for ^-plane 
and /3-plane models, respectively. The position of the center of vortex, obtained nu- 
merically correlates precisely with the theoretical results. The computations were 
made by a modified Lax-Wendroff scheme, the method is second order accurate in 
both space and time variables [27], [39], [20]. We apply the scheme to the system in 
the conservation form written compactly as follows: 
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Suppose that the solution domain in 2D is divided into rectangular cells. Let J7™„- = 
U((i — ^)Aa;, (j — |)Ay, nAt) be the value of U in the center of (xi,yj) at the time 
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level t n . Let us denote 



fi,j = + h-hJ+h + /s+ij-i + fi+UH^ 

1 



(«x/)iJ = ^(/ <+ ^-/ < _x ii )", 

1 

Ay 1 

f y i,j = ^(/iJ+i+ZiJ-i)- 
The explicit two-step scheme is given by 

U?+ 1 = U? J + At(5 x F + 8 y G + S)%K 

The computations were performed on a (240 x 240) uniform grid with the space step 
Air = Ay = 0.64 which corresponds to 12.8 km and the time step At = 0.0005 which 
corresponds to 10 sec of the real time. We use the Neumann boundary condition set 
sufficiently far from the vortex domain. Nevertheless, it is possible to use more sophis- 
ticated non-reflecting boundary conditions [16] and introduce an artificial viscosity to 
damp the oscillations [2"B] . 

5. Steady vortices on the sphere. In Sections [2] and [3] we mentioned that 
all equations, e.g. (|2.ip - (|2.3[) , (|2.6[) . p. II) . (|3.2p . can be written in curvilinear 
coordinates, in this case we mean by derivatives the covariant derivatives with respect 
to the metric of a Riemannian manifold. 

For example let us rewrite system (|2.1[) - (|2.2[) avoiding the tensor notation in 
the spherical coordinates, making clear the influence of curvature of the space, as this 
system usually appears in the geophysical textbooks (e.g. [25j ) : 

„ (*> - („ + _2_}y sin* + -U °Z) . A , 
V i?cos0/ Rcos^dXj 



'2^ + [/ sm «/,+ -— „ =/ 2 , 



V ' RcosS R 




d d U d V d „ . , „ , „ , 

where — = — H = — - H =, R is the radius ol the Farth. 

dt dt R cos (j> d\ R d(j> 

Here we use the established notation (A, cf>) for the geographical longitude and 
latitude, A e (— tt, tt), e (— ■§,§) (we mark coordinates by the character tilde to 
avoid confusion with the moving coordinate system, where the latitude and longitude 
are denoted as (A, (f>)). Further, (U, V) and (/i,/2) stand for the components of 
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Figure 5.1. 
is in radians. 



The periodic field of velocity corresponding to <& 



- cos A cos 0. T/ie space scale 



velocity and the vector of exterior forces F, respectively. The Coriolis parameter is 
I = 2 f2 sin0, f2 is again the vertical component of the angular velocity of the Earth. 
The master equation (|3.15p now has the form (in the moving coordinate system) 



<9 2 $ / ,0* x2 



d\d(j) \ d(j> cos 



i_ e* a \ 

3 2 V 5A ; y 



/ i 02$ a 2 $\ 9$ 9$ sin0,9$, 3 _ o 
\ cos 2 </> <9A 2 9</> 2 / <9A d(j) cos 3 <j> d\ 



(5.1) 



Equation (|5.1j) has a solution of the form 

<1>(A, 0) = $(cos 0(a sin A + 6 cos A) + ft, sin 0), where a, b, h are constants. (5-2) 

In particular, one can take 

<&(A,(/>) = 6cosAcos</>, b = const, (5-3) 

to obtain the zero of the respective velocity field at the origin of the coordinate system. 
Thus, this field is 

u(X, 4>) = — bcos Asin0, v(X, 4>) = 6 sin A. (5.4) 

This field plays the similar fundamental role in the case of the spherical geometry 
as the velocity with linear profile in the case of the plane. It is easy to see that 
(u,v) ~ (~~b(f>,b\) at the point A = <p = 0. To obtain a cyclonic (clockwise) vorticity 
at the origin, we can set b < 0. Fig |5.1l presents the field of velocity for b = — 1. The 
velocity of form (|5 ,4[) was used for testing the numerical methods, e.g. [18] . |23j . 

If a 2 + b 2 0, condition (|3 . 1 6|) can be satisfied only for A = IqL, Iq — const. For 
example, the pressure field ttq, corresponding to l5.4[ is 



/, ,\ ^ 6 cos A cos /,,,„, v „ ^ 

7To(A, 4>) — C (o cos A cos + 2lo) , G = const. (5.5) 

2c 

The expansion at the origin has a form 

M\ 4>) = C- b -^dA (A 2 + ^ + o((A 2 + 02)) . (5 . 6) 

ZCq ZCq 
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Equation (|3.8[) can be also solved here; a rather cumbersome solution is expressed 
in elliptic functions Ep{z, k) and Ep{z, v, k) (see [T]): 

7Ti(t, \4>) = G (cos 4> cos A, 6t + cos0sin2A (E F (z, k) — 2E P (z, k, k))) , (5.7) 

with arbitrary differentiable function G. Here 



z 



z 

Ep(z,u,k)= / ^= — 



2 £2 



1 + cos A cos 4> cos A — 1 j— 

k = v = - -, z = — — — Vk. 

1 — cos A cos <p sin A 

It can be checked that for all G(jii,r]2) having bounded partial derivatives for 
T)i = and all 772 we have 



V7Ti(t,A, 



= 0. 

(A,0)=O 



For the sake of simplicity we take tti = 0. Then the trajectory of the steady vortex 
can be obtained from (|3.10|) . In particular, for ft = (non-rotating case) and Iq = 
we get a stationary exact solution (u(A, 4>), ttq(X, (j))). 

Since system (|5.1[) , (|3.6j) has many solutions dependent on the choice of function 
<&, one can construct a great variety of couples of "cyclonic" and " anticyclonic" steady 
vortices on the non-rotating sphere. 

If Q 7^ and/or m ^ 0, then the position of the respective vortex can be found 
from equation (|3.10l) : 

X(t) +2r2sin(X 2 (i))iX(t)) + c V7Ti(t,x) =0. (5.8) 

x=0 

The discrepancy term for A = 20, sin X 2 (0) is 

Q - F = 2Q ( sin(0 + X 2 (t)) - smX 2 (t)) LX(t) + 



' (5.9) 

2fi(sin(0 + X 2 (i)) - sin(X 2 (0))) Lu(x) + c Vtti(^x) -Vtti^.x) , 

L x=o J 

here x = (A, 4>), u = (u, v). For F = it can be considered as small for small t in a 
small neighborhood of the center of the theoretical vortex. 

For a = b = in (|5.2p condition p. 161) is satisfied for A = 2£1 sin <j) L, the respective 
pressure field is 

7r o (A,0) = --h{h + 2Q) sin 20 + C. 

Thus, we get an exact solution with m — for the rotating sphere, moving according 
to (|5.8p . But evidently it is not a vortex solution, it can be better interpreted as a 
zonal flow. 

We do not dwell here on the computations made for the case of the spherical 
geometry. This is a difficult issue, we reserve it for our future research. Difference 
methods that can be useful here are discussed in [3], [3(3], [2"8] . 
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6. Discussion. The vortex motion is intrinsic to the fluid, both compressible and 
incompressible. There is a huge literature, dedicated to the subject (e.g. [T§] . [3"3" ] . [2"T ] . 
[2"2"]1. especial emphases on vortex of rotating fluid and geophysical applications was 
made in [32] and [13j. In particular, the earliest results concerning vortices in the 
compressible fluid were obtained by C.Chree (e.g. [9]), it is remarkable that the author 
always bears in mind the meteorological context. Recently the vortices in compressible 
fluid were extensively studied by Shivamoggi (e.g. [35]). 

It is interesting that a very complicated systems of equations describing the mo- 
tion of compressible fluid possesses global in time solutions of very simple form. For 
example, in the Euclidean space it is a solution with linear profile of velocity, known 
since Kirkhoff. Meanwhile the solution is not physically reasonable since the velocity 
and the respective pressure field rise unboundedly as the space variables go to infinity. 

The present study was inspirited by a well known fact that the tropical cyclones 
(typhoons) are very stable atmospherical structures and moves as a "rigid body" in 
a background field of pressure. There were many attempts to use this fact for a 
description of the trajectory of the vortex. It is possible to apply a beautiful theory 
of point vortex (point singularity) at a sphere [4] to geophysical problems [TS] , [5] , [5] . 
Further, there is a number of paper, where the problem on the typhoon trajectory 
was solved in assumption that the vortex is a square root singularity, e.g. [7], [10] . 
Nevertheless, the assumption on a singular vortex structure contradicts to the data 
of observations |34) , which evidence the linear profile of the velocity field near the eye 
of typhoon. 

We propose a method of computing the position of the center of vortex. The 
technique can be used for any Riemannian manifold, nevertheless our main interest 
is a sphere and a plane as a local flat approximation of the sphere near a point. We 
consider two types of approximations: the /-plane model and the f3 - plane model. 

We solve a problem of finding the trajectory of a steady vortex moving with a 
bearing field of pressure. Namely, we consider a special kind of solution to the system 
describing the motion of compressible fluid on a rotating two-dimensional manifold. 

Our method consists of several steps. 

a) We use the fact that the vertical scale of the atmosphere is small compared 
with the horizontal one. This allows to average the whole atmosphere model over 
the height and to reduce it to two space dimensions. The vertical structure plays 
the crucial role on the phase of formation and maintenance of the stable vortex. The 
averaging procedure hides these vertical processes, this simplifies the model and helps 
to manage with the problem of the trajectory describing. For the sake of simplicity 
we dwell on the barotropic case. 

b) Further, as is known from experiments, near the eye of cyclone the motion 
is axial-symmetric, the divergency vanishes and the tangential component of veloc- 
ity rises linearly (see [34]). We put the origin at the center of the vortex having 
coordinates (Xi{t), A^(i)) (here we follow [7]) and in the new coordinates (t,xx,X2) 
we perform a procedure that can be called an approximate separation of variables. 
Namely, we separate all terms into three parts: the first depends only on the new 
space variables (xx, X2), the second depends only on t, and the third depends both on 
space and time variables however from some reasons can be considered as small. In 
this connection we use a notion of the 5-approximate solution of the reduced system 
of atmospheric motion in a specific sub-domain of the manifold where the system is 
given. 

c) Given a vortex structure we find the stable component of (re-normalized) pres- 



18 



O. ROZANOVA, J.-L.YU AND C.-K. HU 



sure corresponding to the vortex. Further from a linear PDE we find the bearing 
part of the pressure field. Then we get a second order ODE for the vortex center 
position X(t) (in general case it includes the distance from the center as a small pa- 
rameter). In the case of a sphere we find exact stationary solution of the planetary 
scale (non-rotating case) and study an approximate vortex solution in the rotative 
case. 

We investigate the question on the difference between trajectory of non-localized 
vortex on the I - plane where the exact solution can be obtained and the trajectories 
of localized vortex on I- plane and /3 - plane that seems more realistic. We show that 
the difference is reasonably small. 

We do not stop here on the case of the presence of (turbulent) viscosity in the 
model. This situation can be considered in our framework. The difference is in the 
form of the master equation, in other words, the problem is to find a stationary 
solution to the Navier-Stokes equation. The class of such solution is not so rich as in 
the non- viscous case, we refer for a comprehensive review of this question to [38 . 
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